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We present first results on masses of the scalar and tensor glueballs as well as of the torelon from simulations 
of QCD with two light flavours of Wilson fermions. The gauge configurations of extent 16^32 at (3 — 5.6 and 
K = 0.156,0.157 and 0.1575 have been generated as part of the SESAM collaboration programme. The present 
lattice resolutions correspond to — 2.0-2.3 GeV and ratios m^/mp « 0.83, 0.76 and 0.71, respectively. Studies 
on larger lattice volumes and closer to the chiral limit are in progress. 
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1. INTRODUCTION 

Recently, increased attention has been paid on 
lattice glueball computations . It is demand- 
ing to include sea quarks into such studies as two 
new effects are expected to set in: the scalar glue- 
ball becomes unstable and can decay into two 
tt's (at sufficiently light sea quark mass). Also, 
mixing with mesonic / = states might occur. 
In this case, the glueball will just add an addi- 
tional excitation to the singlet state of the L = 1, 
J^'-^ = 0++ meson quartet (nonet in full QCD 
with three light flavours of sea quarks), which 
does not fit into the constituent quark model. In 
fact, one can equally well speak of a meson with 
valence quarks and sea glue as of a glueball with 
valence glue and sea quarks. Only the couplings 
to certain decay channels can discriminate which 
of these two features dominates in either state. 
Depending on the quark masses, a glueball cre- 
ation operator will project onto all these states 
(tttt and mixed meson-glue states). One aim of 
the present study is to look for unexpected or 
dramatic effects around such thresholds. 

In addition to glueballs, we study so called 
torelon states, i.e. flux tubes wrapping around 
the lattice (with periodic spatial boundary con- 
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ditions). In case of the valence quark approxi- 
mation, the mass of such a state is (up to finite 
size corrections) expected to equal Lsan where 
Ls denotes the spatial extent of the lattice (here: 
Ls = 16) and n — Kja} stands for the string 
tension. When dynamical sea quark flavours are 
switched on, the lightest such state {A\g) can, at 
sufficiently large lattice extent, decay into either 
a scalar meson or into two tt's or mix with glue- 
ball states through intermediate mesonic states. 
In the torelon we expect to be sensitive to "string 
breaking" effects. 

It is well known from quenched studies that 
gluonic observables suffer from large statistical 
fluctuations and small signals Jl|-|^. Early stud- 
ies with staggered sea quarks |^ indicate that it 
might be extremely difficult to obtain significant 
results on glueballs with dynamical fermions at 
all in realistic time. In view of the interesting 
physics involved, we wish to explore the situation 
with our medium size sample of configurations. 

2. METHOD 

We compute matrices of zero-momentum 
projected plaquette-plaquette and Wilsonline- 
Wilsonline correlators at various fuzzing levels . 
Two different fuzzing procedures have been ap- 
plied for this purpose: 

1. successive factor two blocking O] with the 
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straight connection weighted by a coeffi- 
cient a = 1, in respect to the neighbouring 
four spatial staples, 

2. alternating between conventional APE 
smearing steps with a = 4 (keeping the 
length of the fuzzed link constant) and the 
above factor two blocking steps. 

Within procedure 1, the fuzzing levels are labeled 
by n = 1, . . . , 3, corresponding to links of effective 
length 2" in lattice units. Within procedure 2, 
two successive levels correspond to fuzzed links 
of equal length, i.e. levels and 1 yield links of 
extent 1 while 6 and 7 yield links of length 8. 

Procedure 1 has been applied onto 2196 suc- 
cessive trajectories at k = 0.157, V = 16'^32, 
(3 = 5.6 while procedure 2 has so far been ap- 
plied to 204, 274 and 260 configurations, sepa- 
rated by about 10 trajectories, at k = 0.156, 
K = 0.157 and k = 0.1575, respectively. Prior 
to statistical analysis the data has been binned 
into sufficiently large blocks to achieve stability 
of the error estimates in respect to the block size, 

1. e. to avoid autocorrelations. Statistical errors 
have been computed from the scatter between 500 
bootstrap samples in either case. 

Operators with quantum numbers Aig (the 
lightest torelon state), T^'^ (projecting onto the 
jPC ^ 0++ glueball) and E++ (2++) have been 
constructed at each fuzzing level. After diagonal- 
izing the 2x2 submatrix between the two basis 
states with best ground state overlap in each case, 
masses of the glueball and torelon states are ex- 
tracted. 

Applying fuzzing procedure 1, levels 2 and 3 
have been found to project best onto the glueball 
ground states while just level 3 resulted in a sat- 
sifying overlaps for the torelon case. For fuzzing 

2, levels 6 and 5 (7 and 6 for the torelon) are opti- 
mal. After diagonalizing the corresponding corre- 
lation submatrices, we obtain ground state over- 
laps of 72 ± 7% for the 0++ gluebaU and 70 ± 9% 
for the torelon by use of fuzzing 1. With proce- 
dure 2, glueball overlaps of 71 ± 6%, 81 ± 5% and 
79 ± 5% and torelon overlaps of 93 ± 8%, 98 ± 6% 
and 89 ± 4% have been achieved for the three k 
values, respectively, resulting in an earlier onset 
of plateaus in local masses and, thus, reduced sta- 



tistical errors of final results. 
3. RESULTS 





1.4 r 




1.3 - 




1.2 - 




1.1 - 


O 


1 - 


;t-i) 


0.9 - 


o 


0.8 - 


a> 




o 


0.7 - 




0.6 - 




0.5 - 




0.4 - 



Torelon 



16 K 



Figure 1. Torelon effective mass (fuzzing proce- 
dure 1). The rightmost point is the expectation 
from the static potential. 
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Figure 2. Torelon effective mass (fuzzing proce- 
dure 2). 



In Figs. |i| and ^ we compare our results for 
the torelon mass at k = 0.157 from the two 
fuzzing procedures. The numerical values are 
M{Aig) = 0.698^^^ and 0.679l^2> in agreement 
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with an expectation of 16A' — 0.712_j^7, as com- 
puted from the Hnear slope of the interquark po- 
tential 1^. Note, that the leading order finite size 
effect, 7r/(3i|) « 0.004, is small. 



comparison, the quenched data points have been 
scaled by the ratio of the corresponding rp-values: 
5.35/5.55 § within the fi gures. Again, the two 
fuzzing procedures yield consistent values. 
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Table 1 

Torclon masses. 



p = 6.0 



Figure 3. Effective glueball masses (fuzzing pro- 
cedure 1). The rightmost data points are expecta- 
tions from quenched simulations at /3 = 6.0 [3,6]. 
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Figure 4. Effective glueball masses (fuzzing pro- 
cedure 2). 



In Figs. H and ^ we relate our n = 0.157 re- 
sults for glueball masses to those expected from 
quenched simulations at /? = 6.0 i.e. at 

similar lattice spacing (a^^ « 2.1 GcV). For 
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In Table g, we have collected results of the 
torelon mass as extracted by use of fuzzing pro- 
cedure 2 for all three k values. For comparison, 
we have included the expectation from the (ef- 
fective) slope of the interquark potential at large 
distance, k ||^, as well as the mass of two (non- 
interacting) tt's 0. Note, that at k = 0.1575 we 
are traversing the tttt threshold in the sense that 
the TT turns lighter than half the torelon mass ex- 
pectation from the potential. Thus, it is not clear 
whether our ground state still corresponds to the 
torelon or to a two pion state. This might ex- 
plain why the torelon seems to become somewhat 
lighter than the expectation in this case. 



Table 2 

Glueball masses. 
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0.1560 
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1.14(12) 
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In Table || results for the scalar and tensor glue- 
ball masses are displayed. Note, that the value 
m(0++)a = 0.56(6), obtained at k = 0.1575, 
is consistent with 2mT^a = 0.58(1) 0, opening 
this decay channel. In Table |3[ we compare un- 
quenched with quenched results |^ obtained at 

= 6.0 as well as quenched results from Refs. |l|- 
HH] that have been extrapolated to the continuum 
{P — oo) quadratically in a. The scale r^^ from 
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Table 3 

Comparison with quenched resuhs. 
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0.1560 
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3.69(23) 


1.67(15) 
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4.22(14) 


1.40(15) 
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the interquark potential corresponds to 400 MeV. 
All results are in agreement with the quenched 
data at (3 — 6.0. Deviations from the continuum 
limit are expected to depend linearly on a with 
dynamical Wilson quarks. 

4. CONCLUSIONS 

We have studied masses of the scalar and ten- 
sor glueballs as well as of the torelon with two 
flavours of Wilson quarks at (3 = 5.6. All re- 
sults are consistent with those of quenched stud- 
ies at similar lattice spacings. The statistical er- 
rors turn out to be encouraging and will compete 
with their quenched counterparts, once we have 
arrived at final statistics. Between k = 0.157 and 
K = 0.1575, corresponding to ratios mT^jnip = 
0.76 and 0.71 the tt becomes lighter than 
half the torelon mass and string breaking is ex- 
pected to set in on Ls = 16 lattices. Around 
K « 0.1575 the scalar glueball becomes unsta- 
ble, which means that this region deserves fur- 
ther attention. Computations on 24^^40 lattices 
at K = 0.1575 and k — 0.158 are in progress. 
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